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NON-LINEAR  ELLIPTIC  PARTIAL  DIFFERENTIAL  EQUATIONS  AND 
DIFFERENCE  EQUATIONS 

ABSTRACT 

The  Dirichlet  problem  for  the  non-linear  elliptic  partial  differential 
equation 

a(x,y,u(x,y))u,^  +  c(x,y,u(x,y) )u,yy  -  7(x,y,u(x,y) )u  =  0 

is  studied.  It  is  asstuned  that  the  coefficients  are  strictly  positive  and 
Lipschitz  in  the  argument  u(x,y).  It  is  then  proved  that  the  solution  may  be 
uniformly  approximated  by  the  solution  to  the  associated  difference  equation 
provided  that  a  certain  inequality,  relating  bounds  on  the  coefficients,  is 

tottOXSl  XCU* 


3 


This  paper  studies  a  difference  method,  for  a  non-linear  elliptic  partial 
differential  equation  of  the  form 


a(x,y,u(x,y))u,^  +  G(x,y,u(x,y) -  7(x,y/u(x,y) )u  »  0 


(1) 


where  (x,y)  belongs  to  the  open  unit  square  u(x,y)  =  ^(x,y)  for  (x,y)  on 
the  boundary,  of  Q  and  a,  c,  7  are  Lipschitz  functions  of  the  arguments 
x^y^  and  u^x^y) •  The  constant  will  denote  the  maximum  of  the  Lipschitz 
constants  associated  with  a,  c  and  7. 


We  shall  assume 


that  thci e  exist  positive 


UUllO  OCXll  00 


K  Siiiu.  in  Sticli  "uhst 
o 


Kq  >  a(x,y,u(x.y) )  >  m  >  0, 
>  c(x,y,u(x,y))  >  m  >  0, 


and 


Kq  >  7(x,y,u(x,y))  >  m  >  0. 


Place  a  square  grid  on  Q.  The  lattice  points  are  expressible  in  the  form 
(rah,nh)  with  ra  and  n  running  through  a  set  of  integers  and  with  h  denoting  the 
mesh  size.  Let  =  (x,y)  be  a  lattice  point.  The  neighbors  of  P^,  denoted 
by  (v  =  1,  ...  ,4),  are  the  points  of  the  form  (x+h,y),  (x-h,y),  (x,y+h), 
(x,y-h).  Let  consist  of  those  lattice  points  for  which  all  of  their  neigh¬ 
bors  are  in  Q.  Let  3  0^  denote  the  boundary  of  We  place  a  lexicographic 

Order  on  the  points  of  and  denote  these  points  by  P^,  ...  ,  The  points 

of  9  are  denoted  by  Por  a  function  f(x,y),  defined  on 

we  denote  b.y  f.  the  value  of  f  at  P.  .  For  any  M-dimensional  vector  f,  the 
■-  1  . 

i 


or  max 

i;  j 


'"ij' 


I T /~\v»  r\-v  Q  ynoTv^T  V" 

•  V./X  u  mcx  • 


quail  ojr 


M  <.  ! 


UdlV  OCO  L*11C 


In  $1  a  maximum  principle  is  proved  following  the  methods  in  Bers  [!]• 
Existence  and  uniqueness  of  a  solution  of  the  difference  equation  associated 
with  (1)  is  proved  in S 2  under  a  certain  assumption  governing  the  domain  of 
dependence  of  the  coefficients  a,  b  and  7  upon  the  solution  function  u(x,y). 
In  3^  convergence  of  the  solution  of  the  difference  equation  to  the  solution 
of  the  differential  equation  is  proved.  The  condition  imposed  on  the  coeffi- 


CXdJOCj^  uvjii^ci  11X11^  one  uoiuaxii  Ux 


^  *1  <-»  *1  )  I  r  .-P'  v->  •?  j-t 

U.epeiiu.eiiL.e  j  xo  x  eiuwv  etx  xii  uy  u.ex  Xiixilg  a 


new  sequence  of  problems  which  satisfy  the  condition  and  which  converge  to  the 
orginal  problem. 


He  applies  these  methods  to  equations  whose  leading  coefficient  does  not 
depend  on  the  solution. 
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j.  o  ^  IT  j.  rr  ^  j.  ^  f  tt  -  . 

^  ’  ^i^‘'i2  •  -ik'J 
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V-u#-uy 
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§  1.  A  Maximum  Principle 

Let  I  =  (b  )  “"y  M-component  vector =  A  difference  operator  over  fl- 

Ic  ^ 

corresponding  to  (l),  is  ([2],  p.  190) 

T  r.  .  TT  1  _  S  of  a  Mr-  ^IT 

C  ^  a 

where  the  arguments  of  the  coefficients  7j_ 

Let  the  matrix  1(5)  =  (jS^^(I))  defined  such  that  the  (i^k)-th  entry  is 
the  coefficient  of  U,  when  L,  [|.;  U. ]  is  evaluated  at  the  i-th  point  of  the 

i£  il  1  X 

mesh.  The  matrix  -  L(|)  is  monotone  ([3],  P*  ^5)  and  hence  its  determinant, 
det(L(|)),  is  non-zero. 

We  now  state  the  maximum  principle. 

Theorem  1:  Let  U  =  (w  ......  ,U„,U„  =  )  he  any  M-component  vector.  Then 

-  1'  '  N  N+1  M 


1  TT  i 
lu^  lu^l 

17. 

0 

•  •  ^  N  ^  j  “ 

N+1 

T  ft  .  TT  1 


T  * 


f  1  o  ^ 


where  1=1,  . . . ,N,  j  =  N+1,...,M,  and  is  a  constant  independent  of  |  and  U. 


Proof:  Consider  the  case  where,  for  ail  i, 

Lh[li;  u,]  =  0  . 


If 


=  2(a^  +  c^)  +  h  > 


then 


V  -  y  ^  V, 
i  L 


V  =  1 


(1.5) 


where 


A.J  o  —  I.  —  a,  , 

iC,  i<+  i 


1  -  =  >  .  _  =  a. 
■■11  "13  "1 


Since  >  0(  v— 1, . . .  ,^1/  on 


)!  '  !ar!(3^  Q 
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V 

L 

V  =  1 


'"iv 

T 


/  T  7  \  J _ 1  J  —  M 

tiupxxeb  bxiab 

f  <  max  (U. ,  ; . . ,  ,U. ,  )  for  ll.  o 

I  '  ST-'  T  Jl  '  T  _ 

1^  ^  IlllU 


i4^ 

^4' 


1  — 


_ _  /tT  TT  \  TT  ^  r\ 

min  i,  u^.  ^  ,  .  .  .  ,  U^.  ,  J  1  ux-  ^  w 


(1.4) 


The 


equality  in  (1.4)  holds  only  if  for  i  =  1,2,3,  and  7j_U^  =  0, 
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Hence, 


min  (O,  rain  U)  <  (0>  U) 


a  a 


3Q, 


for  i  =  1 , . . . ,  N . 

Consider  the  case, 


i  =  1, . . . ,  N. 


U. 

J 


J  =  N  +  1, . . . ,  M. 


Using  the  definition  of  L  we  have  that 


(1-5) 


(1.6) 


N  M 

-  k.  +  y  (i  =  1,...,  N)  (1.7) 

J  =  N  +  1 

where  a(0  =  (cr^'^(O)  is  defined  in  the  following  way:  a^'^(|)  is  the  coeffi¬ 
cient  of  the  jth  boiondary  point  when  we  consider  L,  [|.;  U.  ].  Note  that  a(0 

n  X  X 

is  rectang\ilar ,  singular  and  the  arguments  are  the  coefficients  a(x,y,|), 
G(x,y,|),  and  7(x,y,^)  evaluated  at  interior  points  of  the  mesh. 


V 

L 

s  =  1 


-IS 


f  t  ^T^ 

V  5  y'-' 


T  T  T 

Letting  U  =  (U^,...,  U^^)  ,  K  =  (kj^_^^,...,  k^^)  and  v  =  (Vj^_^j^, .  . .  ,Vj^)  , 

where  the  superscript  T  denotes  the  transpose,  gives 

L(|)U  =  K  +  a(|)v  . 

Since  L  ^(|)  exists,  we  have 


“  (1.8) 

where  G(|)  =  L"'^(|). 


U  =  L“^(0  K 


1  (w<j(Ov 


or ,  for  i  =  i , . . . ,  N , 


N 


TT 
W  , 
1 


V 

L 

s  =  1 


r<iS  /  _ 


M 

r- 

J  +  1 
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If 


N 


i  vU  -  2^  vw  , 

£  =  1 


then 


max  lul  <  C.(t)  max  |L^[|;  U]  |  +  C^(|)^max  ju] 

■  -  a.  “  " 


(1.9) 


where 


r!  e 
“1"=' 


N 

max  ^  |G^^( E ) 1 

i  Z_.  . 

s  =  i 


M 


CjClj^max  ^  ir'‘'’(Ol 


j  =  N  +  1 


In  order  to  prove  that  €2(0  <  1,  it  will  be  shown  thatp  (|)  and 


i  i 

^  il®  interval  [0,  1]  for  any 

J  =  V+  l' 

Let  U  be  a  solution  of  the  equation 


T  ft  .  rr  1  =  n 


i  =  1, . . . ;  N. 


wiun 


V  =  5  fKronecker  6)  .i  =  N  +  l . M. 

j  j£  '■ 


Then  =p^'^(|)  and  by  (I.5) 


0  <  U.  <1 

—  1  - 


i  =  1,...,  N. 


A 


Leo  U  be  a  soluoiun  of  "the  equation 


L,  [F .  :  U.  1  =  0 
'i-' 


i  =  i  _ .  .  .  w 

—  ^  f  .  .  .  j 


l.TT  +V. 


V.  =  1 
J 


j  =  N  +  1, M. 


Then 


M 


,,  V~'  _i.i/  > 

=  X 

j  =  N  +  1 


and  by  (1-5) 


0  <  <  1 


We 


will  now  obtain  a  boiand  for  C^(|). 


Let 


where  a  = 
Then 


V(x,y)  =  e^ 

-lll/2 


(K^  +  IV 


and  n  is  assumed  to  lie  in  the  strip  0  <  x  <  1. 


0  <  V  <  e 


a 


and 


Lh[5i;  VJ  >  e'^[2a.a^  -  7^1  >  1- 


i  s 

Observing  that  G  (^)  <  0.  since  -L(|)  is  monotone^  (l=8)  gives 


N 

I 

s  =  1 


nlS/. N 


T  r  k  Tr  T  - 

i^h^s  ^  ''s-'  ^ 


M 


j  =  N  +  1 
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Hence,  for  all  ^ 


s 


majc  V 

n 

min  L^[|j  V] 


a 


<  e 


32. 


Existence  and  Uniqueness 

It  has  been  shown  that  any  solution  U  of  the  equation 


J  U.  ]  =  0 
h”"  1  1 


/  O  T  \ 


must  satisfy  the  ineqixality 

max  I U . I  <  max 

i  ^  “  J 

where  i  =  1, . . . ,  N,  v .  =  ^(P.),  and  j  =  N  +  1,...,  M.  By  [1]  there  exists  a 
unique  solution  U  to  (2.1)  for  any  given  4.  Hence  to  every  M- vector  |  there 
corresponds  a  unique  solution  U.  The  dependence  of  U  on  |  will  be  denoted  by 


(2.2) 


u  =  <D(e)  . 


(2.5) 


The  relation  $  is  a  function  satisfying  (2.2).  By  (I.8)  the  function  $  is 

continuous  and,  by  Theorem  1,  can  be  restricted  to  the  N-dimensional  sphere  of 

radius  max  1 v . 1 . 

J  J 

We  are  seeking  a  fixed-point  solution  of  (2.l)  i.e.  solutions  U  for  which 
U  =  o(U).  Application  of  Brouwer's  Fixed-point  Theorem  ([4],  p.  557)  assinres 
the  existence  of  at  least  one  such  fixed-point  of  (2.5). 


Before  proceeding  with  a  uniqueness  proof  the  following  assumption  con¬ 
cerning  our  differential  equation  will  be  made:  Condition  A.  There  exists  a 
set  n'  such  that  W  and  its  closure  are  proper  subsets  of  0  and  in  n  -  O'  the 
coefficients  a,  b  and  7  do  not  depend  on  the  function  u. 


Thpn-rem  2:  Let  a(x,y,u),  b(x,y,u)  and  7(x,y,u)  be  Lipschltz  functions  of  their 


L  WX 


arguments  (x,y,u) 
where 


Let  h  be  so  small  that  satisfies  yGK^m.  ^  <  (5N^C^) 


fv  V  1 

moA 


^  ‘O 

at  most  one  fixed-point  solution 


and  C  =  '^ax  f  1 

0  «.  I  1  - 


1.  Then  there  exists 
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Proof:  Astjurae  there  are  two  fixed-point  solutions,  say  V  and  W.  Then  (1.8) 
gives 

|V.  -  W.|  =  l5,(V)  -  (1>(W)1 

M 

+  1 


V 

L 

^  =  1 


v.yu  -  u  V "  y ;  v  .  I 


M  N 

<  I  |a|  I  -  I  |vl  I  V  VI  g'-^(V)  -  g'-Vw)  I  , 

J - V+ 1  i^l 

where  the  last  inequality  follows  in  virtue  of  Condition  A. 

Since  G  is  the  inverse  of  the  matrix  L,  one  has 

G^^(V)  =  L^^(V)  •  [det[L(V)]]“^ 

where  L^^(V)  is  the  (s,i)  -  cofactor  of  L(V). 

Then 


G""(V)  -  G  (W)l  = 


uc  j-i 


det[L(V)] 

•  det[L(w) 


Sl,„^ 
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Similarly, 


(2.7) 


lj-^J-rT/tr\T  JH  ^4-  r  T  /t.t\  1  1  ^  1  i  T.T  I  i 

jatJULi^V  v;  j  -  ucv.lj_,vw;  J  |  5  I  ''  "  "II  ' 


Substituting  (2.5),  (2.6)  and  (2.7)  into  (2.1^)  yields 


jV.  -  Wj  <  (M  -  N  -  l)N^C^(50K2m  ^)"i  jv  -  Wj 


.0  M. 


If  N  is  chosen  s\iff iciently  large,  then  the  coefficient  of  |  |V 
remains  less  than  one  and  we  have  the  contradiction 


W 


(2.8) 
in  (2.8) 


W 


V  -  W 


Remark:  This  proof  shows  that  the  mapping  (j)  is  a  Lipschitz  function  with 
Lipschitz  constant  less  than  1. 

§  5.  Convergence 

Theorem  3:  If  the  hypotheses  of  Theorem  2  are  satisfied,  the  fixed-point 
solutiion  of  t.h6  differoncG  scLUSt-ion  (1*1)  conY6i*g6s  unifornily  “to  “tho  solution 
of  the  partial  differential  equation  (l). 

Proof;  Assume  the  solution  u  is  completely  known.  Let  be  the  solution 

of  the  difference  equation 


with 


a,D,  +  c.B^  =  0 

i  h,xx  1  1  h,yy  i  i  i 


= 


a  =  a^x  .V  .  )  , 

■'i  "  >  ’ 

/  .,(p-l)  > 

7,.  =  7;,x,  ,y. ,  )  , 


11  1 


n  ^ 


W'.^'=  u.  ,  j  =  N  +1,....,  M,  and 
J  J 

u.  ; 


(5.1) 


here  D. 


V,  U. 

h,xx  1 


etc.  are  used  to  denote  the  coefficients  of 


etc.  in  (l.l). 
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From  (3*l)  and-  (l-O)  we  have 


luff* 

'  1 


1  ' 


y  pij  V. 

" '  (p-2) 


V  . 

1 


(3.2) 


ly(p-l)  _  y(p-2)- 


.4.2 


-2,N 


where  G,  =  5K  C  (50K,m~  )  and  P  corresponds  to  the  n-th  equation  in 

^  o  p  o  1  ^ 

\  —  / 

(3.1).  From  [1]  we  conclude  that 


lw(p)  _ 


^  P?"^  .  £  /'h'l 
^  tv**/ 


f  -z  z\ 


where  €(h)  -» 0  as  h  -»o.  Hence,  the  sequence  J  1  is  a  Cauchy  sequence  and 

('h')  '■  ^ 

has  a  limit  which  we  call  ' .  This  limit  is  a  fixed-point 


soiuuion  oi  our 


difference  equation. 

Repeated  application  of  (3*3)  yields 

llW^P)  -  <Cj^(p)  •  £(h) 


(3.4) 


where  Cj^(p)  =  \  Cy  By  hypothesis  we  have  that  Cj^(p)  is  uniformly 

J  =  1 

bounded  in  p  and  consequently  we  obtain  the  sought  inequality 
,  ,  /'h'l 


U|  I  <  C  .  €(h) 


(3.5) 


where  C  =  lim  C. (p), 

p-»  CO  ^ 


§  4.  Removal  of  Condition  A 

In  |2  a  condition  governing  the  domain  of  non-linearity  of  the  coefficients 
was  introduced.  Here  it  shall  be  shown  how  this  condition  may  be  removed. 

Let  fi  =  {(l  -  e)x,(l  -  e)y:  (x,y)  efi}  with  c  >  0.  Then  and  its  closure 
are  proper  subsets  of  n.  Let  ^  be  extended,  in  a  continuously  differentiable 
way,  into  the  exterior  of  R. 
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Let  a  new  problem  be  posed  as  follows: 

Find  a  solution,  over  Cl,  of  the  equation 

a^(x,y,u)u^^  +  c'^(x,y,u)u  -  7‘^(x,y,u)u  =  0  ,  (l+.l) 

with  u  =  ^(x,y)  for  (x,y)  can.  Here 


))  for  (x,y)  ecf 

X  -  t 

^  1  g))  IX, y;  eS2  -  U' 

and  c*",  7*"  are  defined  analogously. 

Clearly  a^ ,  c",  and  7*  are  Lipschitz  functions  of  their  arguments  (x,y,u) 
and  as  e  0  these  coefficients  converge  uniformly  to  a(x,y,u),  c(x,y,u)  and 
7(x,y,u),  respectively.  Hence  the  solutions  to  (4.1)  will  converge  uniformly 
to  the  solution  of  (l.l)  as  c  0  ([6]  p.  152  and  [4]  p.  597). 

From  §3  there  is  a  uniform  approximation,  U®,  to  the  differential  equation 
corresponding  to  u^(x,y)  i.e. 


Tr  11^  —  ^ 


✓ 

f  a(x,y,u(- 


1  /  //X 

--7  ^ 


ju^  -  u^  I  <  k(c>h) 

where  h  is  the  mesh  size,  e  >  h,  k(e,h)  -»  0  as  h  0  and  IT  is  the  solution 
of  the  e -difference  equation.  Also,  for  the  solution  u  of  (l.l), 

\u^  -  uj  <  £(e) 


where  £(e)  -»0.  Hence, 

|U^  -  u]  <  k(€,h)  +  iiz)  . 

The  following  is  therefore  proved: 

Theorem  4:  If  the  hypotheses  of  Theorem  2  are  satisfied,  if  €  >  h  and  e 
then  there  exists  a  solution  of  the  difference  equation  wMch  converges 
uniformlv  to  the  solution  of  the  differential  equation. 


GREGOR^’  TV  McALLISTER 
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a(x,y,u(x,y))u,^  +  c(x,y,u(x,y) )u,^  -  y(x,y,u(x,y))u  =  0 

is  studied.  It  is  assvuned  that  the  coefficients  are  strictly  positive  and 
Lipschitz  in  the  arigument  u(x,y).  It  Is  then  proved  that  the  solution  may  be  uni¬ 
formly  approximated  by  the  solution  to  the  associated  difference  equation  provided 
i that  a  certain  inequality,  relating  bounds  on  the  coefficients,  is  satisfied. 


a(x,y,u{x,y))u,^^  +  c(x,y,u(x,y) )u,y.y  -  7(x,,y,u{x,y))u  =  0 

is  studied.  It  is  assumed  that  the  coefficients  are  strictly  positive  and 
Lipschitz  in  the  argument  u(x,y).  It  is  then  proved  that  the  solution  may  be  uni- ^ 
fonALy  eipproxijnated  by  the  solution  to  tPie  associated  difference  equation  provided 
that  a  certain  inequality,  relating  bounds  on  the  coefficients,  is  satisfied. 
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is  studied.  It  is  assumed  that  the  coefficients  are  strictly  positive  and  is  studied.  It  is  assumed  that  the  coefficients  are  strictly  positive  and  j 
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is  studied.  It  is  assumed  that  the  coefficients  are  strictly  positive  and 
Lipschitz,  in  the  argument  u(3c,y).  It  is  then  proved  that  the  solution  may  be  uni-  • 
formly  approxiimated  by  the  solution  to  the  associated  difference  equation  provided  ! 
that  a  certain  inequality,  relating  bounds  on  the  coefficients,  is  satisfied.  I 


UNCIASSIFIED 

Dirichlet  problem 
Non-linear  differen¬ 
tial  equation 
E.lliptic  partial  dif¬ 
ferential  equations 
Difference  equations 


